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Abstract. We consider the integrable one-dimensional (^-function interacting Bose 
gas in a hard wall box which is exactly solved via the coordinate Bethe Ansatz. The 
ground state energy, including the surface energy, is derived from the Lieb-Linigcr 
type integral equations. The leading and correction terms are obtained in the weak 
coupling and strong coupling regimes from both the discrete Bethe equations and 
the integral equations. This allows the investigation of both finite-size and boundary 
effects in the integrable model. We also study the Luttingcr liquid behaviour by 
calculating Luttingcr parameters and correlations. The hard wall boundary conditions 
are seen to have a strong effect on the ground state energy and phase correlations in 
the weak coupling regime. Enhancement of the local two-body correlations is shown 
by application of the Hcllmann-Feynman theorem. 
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1. Introduction 

Quantum Bose and Fermi gases of ultracold atoms continue to attract considerable 
interest since the experimental realization of atomic Bose-Einstein condensates (BEC) 
[H El El E] and the pair condensation of fermionic atoms jSl El E|- Particular attention 
has been paid to one-dimensional (ID) Bose gases, which are seen to exhibit the rich 
and novel effects of quantum many-body systems [HI 13 HH EEH EEU 1121 UH US]- As 
a consequence, there has been a revival of interest in the exactly solved ID model 
of interacting bosons. It is well known that the 5-function interacting Bose gas is 
integrable fHl \H\ and can be realized via short-range interactions with an effective 
coupling constant giD This constant is determined through an effective ID scattering 
length aiD ~ a 2 ± /a, where a± is the characteristic length along the transverse direction 
and a is the 3D scattering length. The ratio of the average interaction energy to the 
kinetic energy, 7 = "lf 1D , is used to characterize the different physical regimes of the 
ID quantum gas. Here m is the atomic mass and n is the boson number density. In 
the weak coupling regime, i.e., 7 < 1, the wave functions of the bosons are coherent. 
In this regime, the density fluctuations are suppressed and the phase correlations decay 
algebraically at low temperatures. Thus the ID Bose gas can undergo a quasi BEC. 
However, in the opposite limit, i.e., the Tonks- Girardeau limit 7 ^> 1, the bosons behave 
like impenetrable hard core particles, the so called Tonks-Girardeau gas [TS|. In this 
regime the single-particle wave functions become decoherent and the system acquires 
fermionic properties. 

The ID Bose gas is realized experimentally by tightly confining the atomic cloud in 
two (radial) dimensions and weakly confining it along the axial direction. The motion 
along the radial direction is then frozen to zero point oscillations PHI EDI EH E21 EH], 
making the gas effectively one-dimensional. Anisotropic trapping along the radial and 
axial directions can form either a 2D optical lattice or ID tubes. There have been two 
types of ID quantum gases, the lattice Bose gas EH EI] and the continuum Bose 
gas confined in a harmonic potential along the axial direction [22] (see figure From 
a theoretical point of view, the former is usually described by the Bose-Hubbard model 
while the latter is described by the ID interacting Bose gas. The Bose-Hubbard model 
is not integrable, except for a special case [21], which corresponds to two sites. On the 
other hand, the trapping potential along the axial direction breaks the integrability of 
the ID Bose gas. However, the long-wavelength properties of the ID Bose gas and the 
Bose-Hubbard model can be described by a Luttinger liquid, owing to the universality 
of the low energy excitations, i.e., gapless excitations with a linear low-energy excitation 
spectrum and power- law decay in the correlations |25j . 

The exactly solved ID interacting Bose gas has been extensively studied [2HE3EH 
EH HH USE ED E21 1331 El]- The ground state energy and low energy excitations [IE! EE], 
thermodynamic behaviour [2E] , finite-size effects [SZ] , correlation functions [23 EH] and 
Luttinger liquid behaviour [211 EH E] have been investigated via various methods. The 
signature of the ID Bose gas is strongly influenced by the interaction strength and the 



The ID interacting Bose gas in a hard wall box 



3 



(a) 




(b) 




(c) 




( d ) www 



Figure 1. Schematic diagrams of different models used to describe experiments on 
quasi-lD bosons with contact interaction: (a) Bosons on a line with periodic boundary 
conditions (no potential) - exactly solved via Bethe Ansatz ^J]. (b) Bosons confined 
to a hard- wall box - also exactly solved via Bethe Ansatz (this paper) . Under certain 
conditions the hard walls can mimic a trapping potential, (c) Bosons confined by a 
harmonic trap potential - not exactly solved, but closer to experimental conditions, d) 
Lattice version realized by, e.g., experiments with optical lattices. The corresponding 
Bosc-Hubbard type model is in general not integrablc. 



external trapping potential. The effects of spatial inhomogeneity and finite temperature 
are other considerations to be taken into account under experimental conditions. To this 
end, several approximation schemes have been adopted to describe the main features of 
the ID trapped Bose gas. In particular, the local density approximation [TH H21 EH EI] 
is widely used for calculating the density profiles of bosons and fermions in harmonic 
traps. 

Now for a finite number of bosons and finite system size, boundary effects are 
expected to be pronounced at low temperature jlHlEniEI]- Indeed, significantly different 
quantum effects should be exhibited by a finite number of bosons confined in a finite 
hard wall box. For example, in the weak coupling regime, macroscopic states lie on the 
zero point oscillations as if the system undergoes BEC. The density expectation value 
exhibits Friedel oscillations and the correlation decay is slower than in the periodic case, 
due to the enhancement of the density and phase stiffness. The boundary conditions 
also have an effect on the phase correlations near the boundaries. The ground state 
of the ID interacting Bose gas with hard wall boundary conditions has no momentum 
pairing (with a — k for each k) compared to the period case, because of the missing 
translational symmetry. The hard wall boundary conditions have been experimentally 
realized by square potentials with very high barriers jlHj- Most recently, BEC have 
been produced in a novel optical box trap jlH], in which atom numbers are as small 
as 5 x 10 2 . More experiments in this direction can be anticipated [SHI EH]- These are 
our motivations for studying the ground state properties of the ID interacting Bose gas 
confined in a hard wall box. 

The ID interacting Bose gas with hard wall boundaries was solved by Gaudin in the 
early 1970's jH2|. Gaudin calculated the surface energy via the Bethe Ansatz solution in 
the thermodynamic limit. Very recently, this model was studied via Haldane's harmonic 
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liquid theory ^ZlEOj, which describes the long wave-length properties of the ID fluid 
in terms of the density and phase fluctuations. The correlation functions of the Tonks- 
Girardeau gas have also been studied with hard wall boundary conditions |5Hj . 

The paper is organized as follows. In Section El we present the Bethe Ansatz wave 
functions and Bethe equations for the ID interacting Bose gas with hard wall boundary 
conditions. Details of the Bethe Ansatz solution are given in Appendix A. In Section El 
we derive the ground state energy via asymptotic roots of the Bethe equations in the 
strong and weak coupling limits. We derive the surface energy through the continuum 
integral equations in Section HJ In Section El we calculate the ground state energy 
in the strong and weak coupling limits using Wadati's power series expansion method 
|3*T1 EU E5] . A discussion of the connection between the ID Bose gas trapped by an 
harmonic potential and the exactly solved model is given in Sectional The low-energy 
properties are discussed in Section with concluding remarks given in section |H1 



2. The Bethe Ansatz solution 



The ID quantum gas of N bosons with 5-function interaction in a hard wall box of 
length L is described by the Hamiltonian 



h 2 N d 2 

i=l 1 l<i<j<N 

where the hard walls are defined via the boundary conditions 

V(xi = 0,x 2 ,...,x N )=0, V(x u x 2 ,...,x N = L)=0. (2) 

Here gm = h 2 c/m is an effective ID coupling constant with scattering strength c. The 
wavefunction *f> must be totally symmetric in all its arguments, as required for a bosonic 
system. 

For harmonic trapping along the axial direction, the scattering strength is given 
by c = r^T- The trapping potential should be added to the Hamiltonian (jj) as an 
external field. However, harmonic potentials appear to break the integrability of the 
model. Fortunately the integrability of the model is preserved by the hard wall boundary 
conditions [H2]- This provides us with an opportunity to study the signature of the ID 
Bose gas in a hard wall box in an exact fashion. For simplicity, we set h~ = 2m = 1 in 
the following. 

The explicit solution of the model via the coordinate Bethe Ansatz is described in 
Appendix A. The wavefunction is given by 

* {etki} (x 1 ,...,x N ) = ^e 1 ---e^(e 1 A;p 1 ---e^A ; p^)e i ( ei ^ + ---+ e ^^) (3) 

ei,...,ejv P 

where the sum extends over all N\ permutations P and all signs q = ± (see Appendix A). 
The wavefunction is valid in the domain 0<a;i<...<a;Ar<L and can be continued 
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via symmetry in all coordinates Xj. The wavefunction coefficients A(e\k Pxi ■ ■ -e^k PxN ) 
are determined by the Bethe roots, or wave numbers, fcj via 

N-l N 

A(eik P1 ■ ■ ■ e N k PN ) = (-) p ] \{e x k Pi - ejk Pj + ic)(e i A;p i + ejk Pj - ic). (4) 

i=l j>i 

In the above equation (— ) p denotes a (±) sign factor associated with even/odd 
permutations. The wave numbers satisfy the Bethe equations 

e ,^ = _fi (%-^ + ic)fe + ^ + ie) VJ = lj ... iAr . (5) 
7 = i ( k j -kt — ic) {kj + ki — ic) 

The energy eigenvalues usual given by 

iV 

£ = (6) 

Like the periodic boundary condition case |27j . the Bethe roots fcj are known to be 
real for repulsive interactions (c > 0), but they may become complex for attractive 
interactions (c < 0). Here we consider the repulsive regime. Free bosons are recovered 
for c = 0, i.e., k = irn/L,n £ N (see Appendix A). 

3. Asymptotic solutions to the Bethe equations 

In contrast with periodic boundary conditions, the ground state no longer contains ± 
momentum pairs due to the reflection of quasi momenta at the boundaries. As a result 
the total quasi momentum X^=i ^ s no ^ conserved in this case. We first examine the 
asymptotic solutions of the Bethe equations (JHJ) in the strong and weak coupling limits. 

3.1. Tonks- Girardeau regime 

It is well known that in the strong coupling regime, i.e., 7 3> 1, the ID Bose gas 
with repulsive interaction behaves like a gas of weakly interacting fermions • m the 
limiting case c = 00 the exact solution for periodic boundary conditions and harmonic 
trapping has been given for impenetrable bosons [T51 114j . In the hard wall setting the 
fermionic behaviour can be seen from the ground state energy. Define the variables 
Zj = Lkj/N and 7 = Lc/N. Then for 7 ^> 1/N the Bethe equations (0) can be written 
in the asymptotic form 

' (Zj - Z e ) 2 (Zj + Z e ) 2 \ ^ -A (Zj - Z e ) (Zj + Zff) 



=1 v ' ' ' i=\ e'<e 1 ' 



N 



£=1 



2lE ^( f ^ i£ ) + ( £i ± i£ )^ Vi = l,...,iV (7) 



in which the summations exclude £ = j and £' = j. Here we restrict the solutions to 
Zj > 0. The asymptotic Bethe roots 



N V A 7 
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for the ground state energy follow from the condition that the eqns (J7J) be consistent. 
It follows that in this limit the ground state energy per particle is given by 

E ^\. (N + l)m + l) ( l + 2 Mf2l\\ (9) 



N 6L 2X JX y V L c 

We emphasize that these asymptotic solutions are very accurate for the Tonks- 
Girardeau regime. We will compare the ground state energy Q with numerical solutions 
of the continuum integral equation, which is the hard-wall analogue of the Lieb-Liniger 
integral equation, in Section 0] The explicit form for the wave numbers kj also allows an 
in principle calculation of the asymptotic correlation functions directly from the wave 
function (jHJ). Switching back to real physical units, the ground state energy per particle 
(jHJ) can also be written as 

E h 2 n 2 ( , , 1 



Here the bulk energy eo(7) and the surface energy ef(7) are given by 

eo(7) = £ (i + £) 2 (ii) 
^)=y( 1 + ^)"' < 12 > 

A useful quantity is the ID temperature 7\d = ij§^, which is just the ground state 
energy in different units [22] • We plot Xid obtained from ()10|) as a function of the 
interaction strength 7 in figure 121 for a gas of N = 37 bosons confined in boxes of length 
L = 35.25 /im and L = 32.61 /im. These are the same parameters as in Figure 3 of 
Ref (22] • The dashed horizontal lines are the corresponding values of 7\d in the Tonks- 
Girardeau limit. We remark that for hard wall boundary conditions, the particle density 
distribution is rather flat and homogeneous. It can be seen that Tid increases rapidly as 
7 increases in the weak coupling regime. It then slowly approaches the Tonks- Girardeau 
energy as 7 tends to infinity. In an actual experiment, the length of the atomic cloud 
varies with the interaction strength. In that case the Tid will increase smoothly as 7 
increases. Most significantly, is sensitive to the length of the hard wall box. The 
smaller the length of the box, the larger the 'quasi- momentum' of the particles. We 
note also that the surface energy is positive. In the thermodynamic limit, the ground 
state energy for periodic boundary conditions is proportional to the linear density n. 
Therefore, in the thermodynamic limit, the ground state energy per particle of the Bose 
gas in a hard wall box can be considered as an excited state of a Bose gas with 2N 
particles in a periodic box of length 2L jS2]- We will study the ground state properties 
for the hard wall box further in Section 0] 



3.2. Weak coupling regime 



The ground state properties in the weak coupling limit are both subtle and interesting. 
So far it has proved difficult to reach the weak coupling regime via anisotropic trapping 




Figure 2. The ID temperature Tid obtained from the asymptotic result l|10|) versus 
interaction strength 7 for N = 37 87 Rb atoms confined in ID boxes of length 
L = 35.25/^m and L = 32.61/jm. The asymptotic result (|10|) is valid for a wide 
range of interaction, i.e. 7 3> 1/N. The horizontal dashed lines arc the corresponding 
temperatures in the Tonks-Girardeau limit. The inset shows the ID temperatures in 
the regime 1 < 7 < 30. 



in experiments [22] and it is necessary to sharply define the criterion for weak coupling 
|41| I43j. In the experiment, the ID regime is reached if the radial zero point oscillation 
length l = \JTij (mu±) is much smaller than the axial correlation length l c = h/ y/rnjl, 
where u± is the frequency of the radial trap and fi is the chemical potential of the ID 
system, thus the condition is /i <C ftw± for the ID trapped system. In general, 7 < 1 
at zero temperature is referred to as the weak coupling regime. The Thomas-Fermi 
regime is usually reached for /1 ^> Huj, where u is the axial oscillation frequency for 
the harmonic trap. In this regime the kinetic energy term can be neglected and the 
system has a parabolic density distribution profile [H], referred to as the Thomas- Fermi 
BEC. However, in the regime fi tko the system is considered to have a macroscopic 
occupation in the ground state of the trap with a Gaussian density profile |41j . 

For the hard wall boundary conditions, the leading terms in the ground state energy 
can be obtained through asymptotic solutions of the Bethe equations (jHJ) in analogy 
with the periodic case [Bl]. Here the wave numbers kj for the ground state satisfy the 
algebraic equations 
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Algebraic equations for periodic boundary conditions have arisen in a number of different 
contexts [33], most notably in the integrable BCS pairing models [56J. The roots of 
such equations also describe the equilibrium positions of potentials in Calogero systems 
associated with Lie algebras [57j . 

The solutions of (fTT?|) give the ground state energy 

E n 2 3{N-l)c 

— « h — . (14) 

N L 2 2 L 1 ' 

This is quite different from the periodic boundary case, for which the leading term in 

the ground state energy per particle is -jl ~ (N — l)c/L. Here, due to the reflections 

at the boundaries, the ground state energy for hard wall boundary conditions is larger 

than the energy for periodic boundary conditions. For the integrable ID Bose gas with 

periodic boundary conditions, the ground state energy per particle is known to be given 

by § « §f e (7) with e (7) ~ 7 (l - 5^^) 021 EHl EH • We argue that this result 

holds in the regime 1/N 2 < 7 < 1. The leading term of the ground state energy per 

particle for periodic boundary conditions, i.e. ~ ^~7 holds in the mean-field regime, 

i.e., 7 <C 1/N 2 , for which the correction term is proportional to 7 2 rather than 7 3 / 2 . 

This discrepancy is not totally unexpected, as the Lieb-Liniger integral equation is only 

valid up to terms of order 1/L. If the interaction energy is much smaller than the scale 

of 1/L, i.e., if 7 1/N 2 , results derived from the Lieb-Liniger integral equation are 

no longer valid in this very weak coupling regime. On the other hand, in the regime 

7 > 1/N 2 the zero point oscillation kinetic energy is much smaller than the interaction 

energy and is thus negligible - it is here that one can derive the ground state energy 

asymptotically from the continuum integral equation. Finite-size discrepancies between 

the discrete and integral equation approaches have also been noted in Ref. 



4. The surface energy 

Taking the logarithm on both sides of the Bethe equations (0) gives 
7 T f kj — ki kj + he ) 

kjL = nnij — >^ \ arctan h arctan > (15) 

i=x ^ c c ) 

where j = 1, . . ., N and mj are ordered positive integers, i.e., 1 < m\ < . . . < m^. 

Here for later convenience we have denoted the Bethe roots by k. Our calculation takes 

advantage of the fact that in the thermodynamic limit, the ground state energy of the 

iV boson system in a box of length L is equivalent to one half the energy of 2N bosons 

with length 2L and periodic boundary conditions, as pointed out by |22J. It is thus 

convenient to derive the surface energy from a periodic system of 2N bosons with a 

length 2L, for which the Bethe equations are [To] f 

kjL = 2nlj — j arctan — — + arctan J — - | (16) 



% Of course, one may also obtain the same free energy by treating the Bethe equations Ijl5(l directly, 
see, e.g., Ref. 
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for kj > and Ij are half-odd integers. We now introduce the notation k_j = —kj and 
I_j = —Ij. The difference between kj and kj can thus be written as 

(kj — kj)L = nej — | arctan — — arctan — -1 (17) 

e=-N ^ c c J 

where j = —N, . . . , N and 6j is a sign factor. 
The surface energy is given by 

N 



2 

j=-N 

Using kj — kj < ir/L and taking the Taylor expansion of equations (|17|) at kj = kj + Akj 
gives 



It follows that 



+ (A;,- - kg) 2 ' 
i=-n y ] ' 



Ak< 1 



A c \ _ 1 / ^ cAki \ 

L^^ + ikj-hf) = L [* e * + ^ c* + (kj - W 1 ' 



3 



-N y J w / \ i=-N 



Let us define kj + \ — kj = 2L ^ k ; , where /(&;) is the distribution function ^B], then the 
Bethe equations (fTBj) become 

2,/ W = l + 2c/ B B c2+ g,. )2 ^'- PD 

Here we use the density n = W and define the cut-off momentum 5. Subsequently, 
equation (fT9"j) becomes 

Here e(k) = sgn(fc). Further defining f{(k) = LAkf(k), the surface energy is given by 
Ef= f kf { (k)dk (23) 



J-B 

where ft(k) satisfies the integral equation 

• B /,(*' 



m = ^ k) + iJ^-J&l_*. (24) 

After the same rescaling as introduced in Ref. [TB], i.e., fc = Bx, c = BX, f(Bx) = 
g(x), we find the ground state energy per particle to be of the form (fTUJ) . The bulk and 
surface energies are given by 



7 3 f 1 

e o(l) = / 9o(x)x 2 dx 
ef(7) = y 2 J x 



(25) 



gi{x)xdx (26) 
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where 

/ \ 1 , x f 1 9o(y) , / 07 x 

9 ° {x) = ^ + *L \* + (x- y y dv (27) 

M-fw+y ^+t-v)* * (28) 

with the cut-off condition 

7 J g (x)dx = A. (29) 

There are various methods which can be used to solve the above equations (J27|) and 
()28[) . In the next section we derive analytic results from these equations in the strong 
and weak coupling limits. 



5. Application of Wadati's power series expansion method 

As remarked in references |52] 155] , the Lieb-Liniger integral equation is closely related to 
the Love equation for the problem of a circular plate condensator jHO]. One can obtain a 
series expansion for the ground state energy from the integral equation [31] . This method 
has also been applied to the Yang- Yang integral equations for the thermodynamics [HI] 
and to the Gaudin integral equation for the attractive 5-function interacting Fermi gas 
|55j . For the Bose gas with hard walls, Gaudin [S2] found the leading surface energy 
term in the weakly interacting limit. In this section, we apply the Wadati method to 
obtain the leading terms in the ground state energy from the integral equations (|2*7j) 
and ([2H]). 



5.1. Tonks- Girardeau regime 



In the strong coupling regime the bulk part of the ground state energy per particle is 
given by [ini EH 03 



e (7) 



7 



(30) 



3 V7 + 2, 

which agrees with the asymptotic result (fTTf . Calculating the first two terms in the 
expansion 

(31) 



gf{x) ^a + a x x 
for the distribution function (1251). we find 



a 
ai 



e(x)( 7 + 2) 
27 



2x 5 



3 7 ( 7 + 2) 2 



e (X 7T 



7(7 + 2)' 

This leads to the surface energy 



7T 

ef ~ — 



7 



2 \7 + 2 



2tt s 



7T 2 ( 7 - 2) 

37(7 + 2) ~ 3 7 2 ( 7 + 2) J ' 



(32) 
(33) 

(34) 
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We see that the constant term in the surface energy (|34j) is the same as in (J12)) . but the 
leading correction term differs. Again this is because the continuum integral equation 
derived from the Bethe equations is only valid for terms of order up to 1/L. 

5.2. Weak coupling regime 

In this regime the leading terms of the distribution function are found to be 



term, as for a large number of particles, the contribution from other terms is negligible. 
The ground state energy per particle in the regime 1/N 2 <C 7 <C 1 is again of the form 
()10|) where the leading bulk and surface energy terms are 



So far we have derived some analytic results for the ground state energy of the 
interacting ID Bose gas with hard wall boundary conditions. One may also perform 
direct numerical calculations using the integral equations (J27)l and (J2EJ), as originally 
done in the bulk ^Bj. In doing this we see that the ground state energy (|l()jl. with |TT| 
and (fH?|) . is consistent with the result obtained from the the integral equations (f2Tj) 
and (J28|) for 7 > 1 , with best agreement found for 7 > 5 (see figure OJ) . A comparison 
between the analytic result and numerical calculation for weak coupling is presented in 
Figure |3J A discrepancy between the numerical and analytic results is evident for weak 
coupling. This implies that the next leading term in the surface energy (f?TTj) is necessary. 
As expected, there is a difference between the finite-size results and the limiting curve 
obtained in the thermodynamic limit. 

6. Comparison with the ID Bose gas trapped by harmonic potentials 

The experimental realization of the Tonks- Girardeau gas trapped in harmonic potentials 
has involved the measurement of momentum distribution profiles f20| l2T] , the ground 
state energy j22] and collective oscillations j2Sj- To model the experiments more closely 
it is desirable to take into account the 'soft' boundaries of the harmonic potential rather 
than the commonly used 'hard' boundaries of a box. Unfortunately the axial trapping 
potential breaks the homogeneity of the integrable model. However, if the density 
varies smoothly in a small interval the systems under consideration can be thought of 
as a uniform Bose gas in each small interval fTJ HI] • This quasiclassical approach 
is called the local density approximation and is used to study the density distribution 
profile in cases where the chemical potential is much larger than the level spacing Huu 
in the ID direction. For the equilibrium state the chemical potential of the system in a 
harmonic trap can be taken to be constant. 




eo(7) ~ 7(1 - ^-v 7 ^) 



(36) 



(37) 
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Figure 3. The ground state energy eo + -kef versus the interaction strength 7 for 
N = 10, 20, 30, 40 particles. For each size there is a comparison between the numerical 
solution of the integral equations (|27[1 and (|28|) and the analytic expression for strong 
coupling I|1U|) . with (|llfl and 112|) . derived from the discrete Bethe equations ©. A 
generally good agreement between the numerical and analytic results is visible. The 
lowest curve is obtained in the thermodynamic limit. 



Applying the local density approximation to the ID Bose gas with periodic 
boundary conditions at zero temperature, we have 

lM{n{z)) + V(z) = /i (38) 

where [i(n(z)) is the local chemical potential of the uniform system and V(z) = ^mu 2 z 2 
is the local trapping potential. We make the Ansatz 

li(n(z)) + V{z) = fio, for^l^i? 

n(z) = 0, for \z\>R 1 ' 



where R is the atomic cloud radius given by R = y ^js- The density profile of the 
system can be obtained by using the normalization condition 

n(z)dz = N. (40) 

With help of the analytical expressions for the ground state energy of the interacting 
Bose gas it is now straightforward to derive the density profiles in the Thomas-Fermi 
and Tonks- Girardeau regimes. 
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Figure 4. Ground state energy eo + -j^et in the weak coupling regime versus the 
interaction strength 7 for JV = 10,20,30,40. For each size there is a comparison 
between the numerical result evaluated from the integral equations (|27J) and Ij28(l and 
the analytic expressions for weak coupling, (|36|) and (|37|l . A slight discrepancy between 
the numerical and analytic results appears for weak coupling. This discrepancy 
becomes small if the particle number is very large. The analytic expressions are 
expected to hold in the region 1/N 2 C 7 « 1. 



For the Thomas- Fermi regime, the energy per particle is given by Eq = ^-n{z)c 



and thus 

n(z) = n° TF (l - -jL- ) (41) 
with central density and Thomas-Fermi radius 

Here c = 2/|aiD | • The average energy per particle is given by 



z 2 



N.L Rtf w uv v " 5 V 4m 



1 f RTF Mn^u, 1 f9N 2 u 2 tfc 2 \ 
E TF n- n(z)E (n(z))dz = - . (43) 



In the Tonks-Girardeau limit, the local chemical potential is fi(n(z)) = ^Y^-n 2 (z) 



and the density distribution is given by n(z) = n^ G ^/l — with 



2Nmcu _ 2HN .... 

ir 2 n V mu 
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The average energy per particle in the Tonks-Girardeau limit is E^g ~ \hNuo. The 
density profile has been studied in the whole regime [TJ. The cloud size expands as the 
interaction strength inceases. In the Tonks-Girareau regime, the interaction-dependent 
radius is given approximately by 

r> r> ( -i 32 -Rtg|oid| \ / A -s 

Indeed this would slow down the increasing of the average energy with increasing 
interaction strength in the weak coupling regime if we consider the length of the hard 
wall box varying with 7 via the relation L = 2R. However, further refinements are 
necessary for finite systems, i.e., for a finite number of confined bosons. 

In the previous sections we have discussed in detail the derivation of the ground state 
energy of the ID interacting Bose gas confined in a hard wall box. This integrable system 
is much easier to treat theoretically than the system with harmonic trapping. The 
experimentally measured ID energy has been compared with theoretical curves obtained 
using the local density approximation [22] . However, the theoretical predictions are not 
convincing for a number of reasons. First it is not clear if the quantity 7 avg presented 
in the Figures of Ref. [22] corresponds to the dimensionless interacting strength 7 in 
the uniform Hamiltonian Secondly, the interaction strength region measured, up 
to 7 av g < 6, may be too small to be sure that the Tonks-Girardeau regime has been 
reached. Our theoretical results indicate that finite-size effects induced from the number 
of particles, the system size and the boundary conditions are not negligible in the weak 
coupling and Tonks-Girardeau regimes. 

There is some similarity between harmonic trapping and hard wall box confinement. 
For 7 = 0, the kinetic zero point oscillation energy is \hw for axial harmonic trapping. 

If we confine the ID Bose gas in a hard wall box of length L = 2R , where Rq = \J~^ 
is the characteristic length of the harmonic oscillator, the kinetic energy per particle 

2 

is jf-tiu) for the hard wall boundary conditions, which is much larger than the kinetic 
zero point oscillation energy, for harmonic trapping. For the Tonks-Girardeau regime, 
if we confine the ID Bose gas in a length L = 2i?TG, the ID energy per particle is 
Y^(| + j^)Nhuj, which is almost the same as the average ID energy E^ G fa ~hNu for 
harmonic trapping. 

The boundary effects are more pronounced in the weak coupling limit. If one takes 
the same zero point kinetic energy and ID ground state energy for the hard wall box 
as that for harmonic trapping, the size of the hard wall box for the 7 = and 7 = 00 

limits should be L = \J^j- and L TG = \J 2h ^ 2 {\ + respectively. Recall that 
we plotted the ID temperature as a function of the interaction strength for the hard 
wall boundary conditions with Ltg = 32.61/im and L = 2/?tg = 35.25/xm in figure El 
These numerical values follow on inserting the physical parameters for 87 Rb atoms into 
the above results, with N = 37, 87 m = 0.1454 x 10~ 24 and u> = 2tx x 27.5. In this way 
figure[2]can be compared with the experimental data in Figures 3A and 4 of Ref [22]- As 
observed in Ref. [22] , the radius R of the atomic cloud for harmonic trapping increases 
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rapidly with the interaction strength in the weak coupling regime, causing the average 
energy to increase rather slowly in comparison with the hard wall case. 



7. Luttinger liquid behaviour 

Many one-dimensional models behave like Tomonaga-Luttinger liquids due to the 
universality of the dispersion relation and correlation behaviour. The low energy 
properties are characterized by power-law decay in the correlation functions with gapless 
excitations. A universal description of the low-energy properties of one-dimensional 
interacting systems has been given in terms of harmonic liquid theory [221 • The ID 
Bose and Fermi gases are included in this theory. The Luttinger liquid behaviour of the 
ID Bose gas has recently been studied with hard wall boundary conditions [17] . It was 
found that the particle density exhibits Friedel oscillations with respect to the distance 
to the boundaries. The phase and density correlations are influenced by the hard wall 
boundary effects. Here we examine this behaviour in the context of the integable model. 



7.1. The Luttinger parameters 



The harmonic liquid approach to the low energy excitations of the ID Bose gas is 
described by the effective Hamiltonian [2H1 EDI HZj 



H. 



eff 



\ hv s 



dx 



— n 2 i 

K 



x) + —[d4){x)) (46) 

'o 

where K = \/vj/v N , with v s = ^v^vj, is the Luttinger liquid parameter. Here vj is 
the phase stiffness, vn is the density stiffness and v s is the sound velocity. For the long- 
wavelength density fluctuations the density p(x) = Po + H{x) has small deviations from 
the ground state density po- The boson field operator is defined as ^(x) = \J p(x)e 



where [p(x),e = S(x — x')e The effective Hamiltonian 



is reduced to 



the quantum hydrodynamic Hamiltonian j50J HTj 

hnv 



tf eff = ^M# f (#(?) 

q>0 



ILK 



(N - N ) 



(47) 



with regard to the particle-hole excitation modes. Here u)(q) = v s q for q p an d the 
wave number is restricted to q > 0. N is the total number of particles in the system and 
No is the number in the ground state, with tf(q) the creation operator of elementary 
excitations. The low energy excitations are well described by the effective Hamiltonian 
(|47jl with Luttinger parameters v s and K. We now study the effect of the hard walls on 
the Luttinger liquid parameters. The density stiffness and sound velocity can be derived 
from the ground state energy via the relations [23 HO] 



vn 



L 



d 2 E 
W 1 




d 2 E 



(48) 
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Figure 5. Plot of the Luttinger liquid parameter K as a function of the interaction 
strength 7 for N — 20 bosons with hard wall boundaries (solid line) and periodic 
boundary conditions (dashed line). 



Although the regime 7 <C 1/N 2 is difficult to achieve in experiments, the boundaries 
nevertheless have a significant effect on the Luttinger behaviour in this regime. Using 
the ground state energy ()14j) . we have 

Here the Fermi velocity vp = hirn/m. We see that the Luttinger liquid parameter 



K ps -^L- y 772 + tends to infinity for 7 — > at a faster rate than for periodic 
boundaries, for which K ps tt j '^fj [IQ] (see Figure EJ). The enhancement of K for hard 
wall boundary conditions leads to a very slow decay of the phase correlation |1UJ I4*7| . 
The momentum distribution exponentially decays as n{p) ~ p _/3 , where (3 = 1 — ^ |47j . 
The system behaves like a BEC in this regime. 

In the Thomas- Fermi regime 1/N 2 <C 7 < 1 we use the ground state energy in 
Section 5.2 to obtain the parameters 



tt 2 V 2vr iV^/7 



vn = v F — 2 ( 1 - + ^= ) , (50) 



^O-i^ 1 ' (51) 
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Figure 6. Plot of the Luttinger liquid parameter K versus the interaction strength 
7 for N = 20 bosons with hard wall boundaries (solid line) and periodic boundary 
conditions (dashed line) in (left) the Thomas-Fermi regime and (right) the strong 
coupling regime. In the Tonks- Girardeau limit K tends to the value 1 + 



A 



7T 



1 



2tt 



In the strong coupling regime we use the result Q to obtain 

-4 



V N = V F 1 + 



1 



1 

2iV 



1 



v F 



1 + 

2A^ 



1 + ? 



2 ' 



A' 



2 

1 + - 

7 



1 + J - 

2JV 



H 


i 


M) 




2N 





(52) 



(53) 



(54) 



(55) 



We show a comparison of the Luttinger parameter A for the different boundary 
conditions in the Thomas- Fermi regime and the strong coupling regime in Figure H3 
The parameter A varies from infinity to 1 + jU. For weak coupling, A increases more 
quickly than in the periodic case. In the strong coupling limit, A tends to 1 for both 
cases. 

The density fluctuations are suppressed in the weak coupling limit, while they are 
enhanced in the strong coupling limit. This can be seen directly from the leading term 



The ID interacting Bose gas in a hard wall box 



18 



of the density expectation value [T7| 

1 



<P(*)> 



n 



1 



7T 



7T 



n2L|sin|f 



sin(27T7T,x) 



(56) 



with a <C x <C L 
coupling we find a 



a, where a is the cut-off length to the boundaries. For weak 
: — ^= while in the strong coupling limit a ~ — . In the weak 
coupling regime, the density fluctuations are suppressed due to the coherence of the 
wave functions. However, in the strong coupling limit, density fluctuations are evident 
due to the decoherence of the wave functions. These effects can be seen directly from 
equation (JHEJ)- 



7.2. Local correlation g2 

It is well known that the local correlations in the ID Bose gas decay algebraically 
[SHI IS]- As for the periodic case, we can calculate the two-body correlation functions 
through the ground state energy for the hard wall boundary conditions. Using the 
Hellmann-Feynman theorem, the 52(7) correlation function is given by [%01 

02(7) = n 2 ^- (0(7) + ^rMl) J • ( 57 ) 



<97 V N 

We thus obtain the correlation function g 2 in the various regions as 



47H n. 



2„2 



[l + ^) , for 7 > 5 



9 2 (7) = {Wl-¥ + ^ , for^« 7 «l (58) 



f™ 2 , for 7 < ^2- 

These results are to be compared with the periodic case, for which 



92(1) = < n 2(i-hay for _i_ « 7 « 1 (59) 



n 2 , for 7 < ^2- 



The enhancement of the correlation function g 2 by the hard walls is largest for weak 
coupling, as can be seen in Figure This is due to backward scattering increasing the 
probability of two particles scattering in comparson with only forward scattering for the 
periodic case. 



8. Conclusion 



We have considered the integrable interacting ID Bose gas in a hard wall box. The 
exact Bethe Ansatz solution for the wavefunctions and eigenspectrum has been outlined 
in Appendix A. The ground state energy, including the bulk and surface energy have 
been derived from the Bethe equations (JHJ) in different regimes. For N bosons, these are 
(i) the mean-field regime 7 <C 1/N 2 , where the interaction energy is much smaller than 
the kinetic energy, (ii) the Thomas- Fermi regime 1/N 2 <C 7 <C 1, and (iii) the strongly 
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Figure 7. The normalized local correlation function <?2 versus coupling strength 7 for 
(left) weak coupling and (right) strong coupling for N = 30 bosons. The enhancement 
of the correlation function by the hard walls is largest for weak coupling. 



interacting To nks- Girardeau regime 7> 1. These results have been compared with the 
ground state energy obtained from the continuum Lieb-Liniger-type integral equations 
in the thermodynamic limit. The latter results, (|26 )) -([28| ) . are in agreement with those 
found by Gaudin |H2]- The emphasis of our approach has been on finite systems and 
the effects of the hard wall boundary conditions. It is seen that the finite-size results 
compare well with those from the continuum integral equation, with the exception of 
the mean-field regime, where the integral equation is not expected to hold. 

A connection to the ID Bose gas trapped by a harmonic potential has also been 
made. The Luttinger liquid parameters, such as the density stiffness, sound velocity, 
and the local correlation function g 2 have been calculated from the ground state energy 
in the various regimes. It is clearly seen that the hard wall boundary conditions have a 
larger influence on the phase correlations in the weak coupling limit. The enhancement 
of the Luttinger liquid parameter K strongly suppresses the fluctuation in the ground 
state density expectation value. The local correlation g 2 is enhanced by the hard wall 
boundary conditions. A significant effect of the hard wall boundary conditions is that 
the wave-like properties of the bosons become more pronounced in the weak coupling 
regime. Significantly, the ID interacting Bose gas confined in a hard wall box can be 
experimentally realized. Future experiments, highlighting the subtle interplay between 
system size and boundary effects in ultracold quantum gases, are eagerly awaited. 
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Appendix A. The Bethe Ansatz solution 

In this Appendix we present the coordinate Bethe Ansatz solution of the ID interacting 
Bose gas with hard wall boundary conditions. The general form of the wavefunction is 
that of the spin-| Heisenberg chain with open boundaries and additional surface terms 

EH- 

First we consider the standard one particle in a box case, N — 1, for which the 
wave function is 

= A(k)e ikx - A(-k)e~ ikx . (A.l) 

From the hard wall boundary conditions (J2J), we have A(k) = A(—k) and e l2kL = 1. 
The energy is E = k 2 . Thus, up to an overall constant, the wave function is 
ty(x) = A(k) sinkx, with quasi momentum k = mr/L, where n is a non-zero integer. 
This describes a trapped particle which has a discrete energy spectrum. If L is infinitely 
large, the particle becomes free and the enery levels will be continuous. 
Next we consider N = 2, with wave function Ansatz 

(zi, x 2 ) = A(kx, k 2 ) e l{klXl+k2X2) + A(k 2 , h) e ^ Xl+klX ^ 

- A{-k u k 2 ) jl-hxi+toxa) _ A(-k 2 , fci) e i t- fcaX1+fclX3 ) 

- A(h, -k 2 ) e i( - klXl ~ k2X ^ - A(k 2 , -h) e ^ k2Xl - klX2) 

+ A(-k t , -k 2 ) e - i ^ +krx ^ + A{-k 2 , -k x ) e - i (^i+^^)(A.2) 

with the domain < Xx < x 2 < L. For x\ ^ x 2 , H^f(xx,x 2 ) = (kf + k 2 )^(xi,x 2 ). For 
%\ — %2i the consistency condition for the wave function to be continous is 
d d \ 

^[xi,x 2 )\ xl=X2 = c*(z 1 ,x 2 ) (A. 3) 



dx 2 dxi 
which leads to the relations 

A(fc 1 ,fc 2 ) = f 1 ~f 2 + iC A(fc 2 ,fc 1 ), 
k\ — k 2 — ic 

A{-k u k 2 ) = A(k 2 ,-h), 
h + k 2 + ic 

mi i \ h + h + ic A . 7 

A(h,-k 2 ) = —— —A(-k 2 , kx), 

ki + k 2 - ic 

A(-k u -k 2 ) = k l 1 ~ ] ?~ 1C A(-k 2 , -kx), (A.4) 

K\ — K 2 + IC 

between the coefficients A(kpx, kp 2 ). Using the hard wall boundary conditions (j2J), the 
unnormalized wave function is 

V(x 1 ,x 2 ) = (k x -k 2 + ic) (h + k 2 - ic) e i{klXl+k2X2) 

+ (kx -k 2 - ic)(kx + k 2 - ic) e i{k2Xl+klX2) 



(A.5) 
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- (fei -k 2 + ic)(fci + k 2 - ic) e ^~ klXl+k2X2) 

- (fei -k 2 + ic)(h + k 2 + ic) e^- k ^ 

- (ki + k 2 + ic)(fci - fez - ic) e 1 ^ 1 * 1 -* 8 * 3 ) 

- (fei - k 2 - ic)(fei + fe 2 - ic) e i (- fca<Bl+fcia,a ) 
+ (fei - fe 2 - ic)(fei + fe 2 + ic) e ~ i(fci:ri+A:2:C2) 
+ (fei - fe 2 + ic)(fci + fe 2 + ic) e" i(fc2:Cl+feia;2) , 

provided that the Bethe equations 

e i 2fcl i = (fei -fe 2 + ic)(fei + fe 2 + ic) 

(fei - fc 2 — ic)(fei + fe 2 - ic) ' 
e i 2fc2 i = (fe 2 -fei + ic)(fe 2 + fei + ic) 

(fc 2 - fei - ic)(fe 2 + fei - ic) ' 
are satisfied. 

Observe that for c = 0, the wave function (jA.5|) reduces to the standing wave 
ty(xi,x 2 ) = sin kiXi sin fe 2 x 2 + sin fe 2 xi sin feix 2 , with fei = niir/L and fe 2 = n 2 ir/L for ri\ 
and n 2 non-zero integers. When c increases, fei and fe 2 also increase as if the boson mass 
increases. The standing wave properties are gradually lost as the interaction becomes 
stronger. 

In a similar way, we can derive the iV-particle wave function given in (JHJ). The 
coefficients are connected to each other via 

, dik>i €-jkj ~\~ IC 



(A.6) 



v4(. . . , 6jfej, 



, 6jkj, 



-A(. 



i tiki) 



.), (A.7) 



for i < j. Application of the boundary conditions leads to the explicit form of the 
coefficients given in (jlj), along with the Bethe equations (j2J). The wave function is 
unnormalized. 
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